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ABSTRACT 



Preliminary studies have been carried out for the 
frequency and intensity stabilization of the output of a 
Sylvania Model 941 Carbon Dioxide laser, by feedback 
control of resonator length. Fluctuations due to discharge 
current and temperature and flow rate of coolant were 
studied. Current regulation of the power supply was found 
to be necessary for intensity stability. Cavity length 
tuning allowed vibration- rotation line selection, with a 
strong hysteresis effect. A computer analysis of error 
signal production for a pyroelectric detector and for a 
HgCdTe detector indicated that phase sensitive detection 
of the first harmonic signal component gave the typical 
frequency discriminant. Operation of a servo loop integrator 
for the feedback loop has been considered. A suitable Stark 
Effect absorption cell for feedback signal generation has 
been constructed, but has not been applied to the laser. 
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I. INTRODUCTION 



The Sylvania Model 941 CO^ laser was purchased to do 
atmospheric propagation studies. It was determined then to 
be too unstable for such studies. A program was undertaken 
to study the instabilities and to stabilize the laser both 
in frequency and intensity. The novel method of Stark cell 
stabilization was chosen as the best route to take. 

Stark cell stabilization, being relatively new, has not 
come into commercial production but has been demonstrated as 
feasible in the laboratory and used on a few prototype 
systems. As a result there is very little information in 
the literature as to the exact details of the technique. 
Particularly, it was required to know how to produce an 
error signal and which detector would be most suitable for 
the stabilized system. 

This thesis describes the preliminary accumulation of 
knowledge required to stabilize the Sylvania laser using 
the Stark cell technique. 
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II. THEORETICAL 



A. ROTATIONAL -VIBRATIONAL SPECTRA 

The quantum states of a molecule arise from three types 
of motion : electronic, vibrational, and rotational. The 

electronic charge cloud surrounding the molecule may be 
represented mathematically by wave functions. This wave 
function representation results in quantized energy states 
just as with a single atom. Motions or changes in the 
electronic charge distribution may result in a changing 
electric dipole moment in which case radiation is either 
absorbed or released. Electronic transitions correspond to 
frequencies primarily in the visible and ultraviolet 
portions of the spectrum. 

The atomic nuclei in the molecule themselves are not 
stationary, but vibrate about some equilibrium position. For 
each electronic level there may be several modes of vibration 
which are also quantized. These vibrational levels result 
in a splitting of the electronic levels. Transitions be- 
tween vibrational levels correspond to frequencies in the 
near and middle infrared portions of the spectrum. Vibra- 
tional transitions which result in a changing electric dipole 
moment occur via absorption or release of photons. 

The molecule as a whole may also rotate about various 
axes. This rotational motion is also quantized and gives 
rise to a splitting of each vibrational level into rotational 
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sublevels. Transitions between rotational levels of the 
same vibrational state correspond to far infrared and 
microwave frequencies. As in the other two types of motion 
only those transitions which result in a changing electric 
dipole moment will occur with absorption or release of 
radiation at these frequencies. The above described energy 
level scheme is represented by Fig. 1. 

Carbon dioxide is a linear, symmetric, triatomic molecule 
It can vibrate in three different modes as pictured in Fig. 2 
These are: v^ , the symmetric stretch mode in which the two 

oxygen atoms move in opposite directions along the longi- 
tudinal line of symmetry while the carbon atom remains 
stationary; v^, the bending mode in which the oxygen atoms 
remain stationary and the carbon atom moves along a line 
perpendicular to the longitudinal line of symmetry; and v_, 
the asymmetric stretch mode in which the oxygen atoms 
remain stationary and the carbon atom moves along the longi- 
tudinal line of symmetry. Only the v ? and v, modes give 
rise to an electric dipole moment and can thus couple 
directly to electromagnetic radiation. The v^ motion 
appears strongly in Raman scattering. 

The v 7 mode can occur equally in two different 
orthogonal directions . Thus the v 7 mode has a twofold 
degeneracy. The quantum number Z is a measure of the 
angular momentum about the longitudinal axis as a result 

of the v^ motion. The vibrational quantum states are then 

Z 

denoted by [v^ v ^ v-] . 
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The vibrational energy associated with a given state is 
approximated to first order by 

E v = ( v i + 7 ) h cv x + (v 2 + j) hcv 2 + (v 3 + j) hcv 3 (1) 

where = 1388.3 cm 1 , \> 2 - 667.3 cm 1 , v 3 - 2349.3 cm 1 , the 
frequencies in wave number units associated with the v^ , v 2 , 
and v 3 modes, respectively. v^ , v 2 , and v 3 in the above 
equation are the vibrational quantum numbers and take on 
values 0 , 1 , 2 , ... . 

The C0 ? energy levels of principal interest in laser 
action and their corresponding frequencies are: [00°1] at 

2349.3 cm 1 , [10°0 , 02°0 ] j at 1388.3 cm 1 , [ 10° 0 , 0 2 0 0 ] j j at 
1285.5 cm 1 , and [01 1 0] at 667.3 cm 1 . The [10 o 0,02 o 0] I and 
[10°0, 02 o 0]jj representations follow the notation used by 
Tyte [Ref. 1] and signify the fact that the [10°0] and 
[02°0] levels are in Fermi resonance. The main transitions 
of interest and these four levels are represented in Fig. 3. 
In the figure the 2 and n notations refer to the electronic- 
vibrational angular momentum. This vibronic angular 
momentum K takes on values 0, 1, 2 ... corresponding to 
5, n, A ... , respectively and is the vector sum of the 
vibrational angular momentum about the internuclear axis 
and the electronic orbital angular momentum about the 
internuclear axis. The 2, n, and A states of a molecule are 
analogous to the s, p, and d states of an atom. The g 



11 



subscript denotes that the vibronic eigenfunction remains 



unchanged upon reflection at the center of symmetry, i.e., 
the eigenfunction is of even symmetry. The u denotes an 
eigenfunction of odd symmetry. The plus superscript denotes 
that the eigenfunction remains unchanged when reflected at 
any plane containing all of the nuclei [Ref. 2]. 

The C0 o molecule can rotate about any axis perpendicular 
to the longitudinal axis. The energy of rotation is given to 
first order by 



where I is the moment of inertia and J is the rotational 
quantum number. From Eq.(2)it is observed that the energy 
difference between two adjacent rotational states is 



where J is the quantum number of the lower state. Substituting 
into (3) the equation 



■fc 2 

Ej = J(J + 1) Ij- J = 0, 1, 2, . 



( 2 ) 




(3) 




( 4 ) 



solving for 1/A the reciprocal wavelength, and taking the 
difference in wave numbers between the spectral lines one 



( 5 ) 



A ,1. _ i i = h at 
A V X 1 ' 2rcT AJ • 

From (5) it is clearly observed that for two adjacent lines, 
i.e., AJ = ±1, the frequency difference is a constant if I is 
a constant. In general, however, I is not a constant, but 
changes slightly with different vibrational states. This 
is the reason for the observed slightly nonuniform spacing 
of the rotational lines. Herzberg [Ref. 2] gives various 
values of the rotational constant for the corresponding 

"fc i 

vibrational levels. A typical value for ^ j— is 0.3899 cm 
for the [ 0 1 1 0 ] vibrational level. This nonuniform spacing 
also shows up in the departure from parabolic of the potential 
curve in Fig. 1. 

The selection rules originate from a calculation of 
the electric dipole moment matrix element taken between the 
initial and final states; 

Pfi = |/^£ erif^dxl . (6) 

Transitions between states which yield the above integral 
zero are forbidden. For a harmonic oscillator the integral 
yields the selection rule Av = ±1. However, as the vibra- 
tional energy increases and/or the molecular structure 
becomes more complex, the vibration departs from simple 
harmonic and Av = ±1 cannot be expected to hold. From (6) 
also comes the selection rule for C0-, that a u state 
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must go to a g state and vice versa. But u to u and g to g 
transitions are forbidden. This is clearly seen in Fig. 3. 

In general the total energy of a molecule to first 
order approximation is the sum of the energies due to the 
electronic, vibrational, and rotational states; 

E ■ E e * E v * E J W 

The I.R. transitions of C0 7 are of the vibration-rotation 
type. This is represented schematically in Fig. 4 for 
transitions between two vibrational levels, v^ and v^. The 
rotational selection rule AJ = ±1 causes the spectrum to be 
divided into two branches: a "P branch" where AJ = +1 and an 

"R branch" where AJ = -1. For CO^ this is a S + - 2 + vibronic 
transition. Thus AJ = 0 is forbidden, which would otherwise 
yield a Q branch at the band center frequency \> Q . The 
frequencies for the P branch are represented by 

v s v JL_ j J = 2, 4, 6 . . . (8) 

0 4tt“ I 

and for the R branch 

v = v + -X- (J + 1) J = 0, 2, 4 ... (9) 

0 4tt l 



where v Q 





TT 



the band center frequency and J is the 
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Table 



rotational quantum number in the lower vibration level. 
Equations (8) and (9) are easily derived from Eq. (2). 

Note the values which the quantum number J may take on. 

CO 2 is a linear symmetrical (point group D^) molecule. 

The two identical oxygen nuclei in the CO^ molecule have 
zero spin and obey Bose statistics. Thus all antisymmetric 
levels have zero statistical weight and the odd rotational 
levels of a Z* state are missing as well as the even 

o 

rotational levels of a Z* state. Consequently the Z* to Z* 
transitions for the 10.4 pm and 9.4 pm bands have the simple 
structure illustrated in Fig. 4 consisting of an R branch 
and a P branch with alternate lines missing. A particular 
vibration-rotation line is denoted by P(J) or R(J) where J 
is the quantum number of the rotational level in the lower 
vibrational state. A few of the more important lines in the 
10.4 pm band along with their wavelengths are given in 
Table 1. Of these lines the P ( 2 0 ) transition is the strongest. 

B. EXCITATION MECHANISM 

The excitation of C0 ? into the [00°1] level is readily 
accomplished through transfer of vibrational energy from 
nitrogen molecules in the v = 1 vibrational level of their 
ground electronic state. From Fig. 3 it is seen that the 
v = 1 excited vibrational level of N^ is about 19 cm'*' below 
the [00° 1] level of C0 2 . The kT average thermal energy of 
molecules at room temperature is about 210 cm'*'. Thus these 
two levels are in almost exact coincidence and a high cross 
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section exists for vibrational energy transfer through two- 
body collisions. This may be represented by the equation 

N*Cv = 1) + C0 2 (00°1) - N 2 (v = 0) + C0 2 (00°1) - 19 cm 1 (10) 
* 

where N 2 represents a vibrationally excited N 2 molecule in 

the ground electronic state. 

* 

N 2 is easily produced by a discharge in low pressure 

nitrogen. The homonuclear molecule N 2 has no permanent 

electric dipole moment; thus av=ltov=0 transition via 

* 

electric dipole radiation is forbidden, causing N 2 (v = 1) 

a 

to be metastable. The downward transition for N 2 ( v = 1) 
is governed almost solely by collisions with other molecules 
and the container walls. 

It should be noted that the excitation cross section 

* 

of C0 2 to the [10°0, 02°0] level by N 2 is very small as the 
difference in energy is more than 900 cm'*'. This fact along 
with the shorter lifetime of the [10°0, 02°0] level allows a 
population inversion in C0 7 to be achieved easily. Helium 
is often added to increase the relaxation rate of the lower 
level . 

C. LINE SHAPES 

Analysis of a spectral line of an atomic or molecular 
transition reveals it to be composed not of a single frequency 
only, but that the line has a finite width or frequency spread 
with an intensity profile which is either Lorentzian or 
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Gaussian in the limit. The frequencies form a continuum from 
zero to infinity with the maximum intensity at v , the 
fundamental frequency of the transition. These two pro- 
files are compared in Fig. 5. 

The Lorentzian or natural lineshape is characteristic 
of homogeneous broadening. In homogeneous broadening, the 
absorbing or emitting molecules are indistinguishable as 
they all absorb or emit at the same frequency. The three 
most common causes of homogeneous broadening are: (1) the 

finite spontaneous lifetime of an excited state, (2) colli- 
sions of molecules with phonons, and (3) collisions of 
molecules with molecules in a high-pressure gas which 
effectively decreases the lifetime. 

The Gaussian lineshape is characteristic of inhomo- 
geneous broadening. Inhomogeneous broadening occurs when 
the absorbing or emitting molecules are distinguishable as 
a result of their absorbing or emitting at different 
frequencies. The two most common causes of inhomogeneous 
broadening are: (1) impurities and variations in a host 

crystal affect the energy levels and thus change the transi- 
tion frequency of nearby emitting molecules and (2) the high 
velocities encountered in hot plasma gases cause the 
frequencies to be Doppler-shifted. Obviously for a C0 ? 
laser pressure broadening will predominate at high gas 
pressures, resulting in a Lorentzian lineshape, and Doppler^ 
broadening will predominate at low gas pressures, resulting 
in a Gaussian lineshape. 
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LINESHAPES 
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FIGURE 



Consider cause No. 1 of homogeneous broadening above. 

It is well known from Fourier transform theory that a pure 
single frequency implies a transition which occurs over an 
infinite amount of time which, of course, does not occur. 
The Lorentzian lineshape may be derived by assuming the 
emitter is a damped harmonic oscillator with an exponential 
decay time t [Ref. 3] . 



c(t) 



E o e 



- t/t 



cos 



2ttv t 
o 



( 11 ) 



The Fourier transform of (11) is equivalent to the convo- 

~ t / T 

lution of the Fourier transform of e and the Fourier 

transform of E Q cos 2TT\> o t: 



F{c (t) } = F{e' t/T } * F{E cos 2ttv t} 
s y o o 



r/°V t/T cos 2irvtdt - i/°°e' t ^ T sin 2iTvtdt] * 

O 0 



E E 

6 (v - v Q ) + j- 5 (v + v Q )_ 



1/t 

— : — t ZT 

(l/x) “ + 4 it v 



i2TTV 

2 7 2 
(1/t)" + 4i"v J 



E E 

<S(v - v q ) + 6(v + v q )_ 



( 12 ) 



where the lower limit of integration is taken as zero rather 
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than negative infinity to correspond to the start of the 
emission. Since negative frequencies are being ignored in 
this analysis, the delta function at v = +v q in (12) may 
be taken with height E q and convolved with the transform of 
e" t/x : 



* [E q <5 (v - v Q ) ] (13) 

Equation (13) yields: 



F{ e (t) } = 



1/t 

(1/t)“ + 4tt 



i2irv 

(1/x) + 4 tt v 



F{e(t)} = 






i 2 it (v - v )E 
v o' O 



(1/t) 2 + 4tt 2 (v - v Q ) 2 



(1/t) 2 + 4tt 2 (v 



v ) 

0 J 



(14) 



Taking the square of the absolute value of (14) (which 
corresponds to the power spectral density in communications 
theory) yields: 



| F{e (t) > | 2 



(1/ t) ~E q 2 + 4 tt 2 ( v - v 0 ) 2 E 0 2 
[(1/t) 2 + 4tt 2 (v - v 0 ) 2 ?“ 



(1/t) 



E 



2 

o 

— 2 

4tt (v - v ) 
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Normalizing (15) according to_^°°g(v)dv = 1 yields 



( 15 ) 
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(16) 



g(v) 



(1/tJ +" 47T ( V - V ) 2 




This is the Lorentzian lineshape in terms of x. Defining 
Av as the full width at half maximum the Lorentzian may be 
found in terms of Av . From the definition it follows that 



Substituting in the proper values of (16) into (17) and 
carrying out the algebra yields 



Equation (18) is worthy of closer examination at this 
point. It is noted that Av is a measure of the frequency 
spread or finite line width. x is the transition time. As 
x decreases the frequency spread Av must increase. Only 
if x were infinite would Av be zero implying a pure single 
frequency transition. This is consistent with the remark 
made earlier concerning transitions occurring in an infinite 
amount of time. 

Substitution of (18) back into (16) yields the Lorentzian 
lineshape in terms of Av : 



g (v Q ) / 2 = g(v Q + fL) . 



(17) 




(18) 
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g(v) 



Av 



(19) 




From (18) and (19) the following is observed to be true: 

Av = — = — l r (20 

ITT Trg(v o ) L 

Equation (20) says that as Av increases or t decreases the 
intensity at v q must decrease. But if the intensity 
decreases at v Q it also decreases at all other frequencies 
within a region near v , frequencies at which a laser will 
normally oscillate. This is one possible explanation for the 
typically observed decrease in power output from a laser 
with an increase in temperature. An increase in temperature 
would increase the pressure without increasing the number 
density and pressure broaden the gain profile resulting in a 
decrease in output at the cavity oscillation frequency. 

Reference 4 gives the time between collisions in a gas 
at pressure p for a molecule of mass M and a diameter a as 




( 21 ) 



The universal gas law is given as 
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p = NkT 



( 22 ) 



where N is the number density of molecules in a gas. 
Substituting (22) into (21) and defining the collision cross 



Substituting (23) into (20) then gives the full width at 
half maximum due to pressure broadening 



Equation (24) clearly shows that Av increases as the square 
root of the temperature. 

Consider now cause no. 2 of inhomogeneous broadening. 
The Gaussian or Doppler broadened lineshape is easily 
derived by beginning with the Maxwell velocity distribution 
function in one dimension which is the probability that a 
molecule of mass M in a gas at average temperature T will 
have an x velocity component v 

X 



? 

section as Q = ira yields 




(23) 




(24) 




r -Mv 2 1 



*(v x ) 



X 



(25) 
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The Doppler-shifted frequency is given as 




( 26 ) 



which upon rearranging gives 



v 

x 




(27) 



Since v is also a random variable and <£(v ) is a probability 
distribution function, a change of variables can be made 
according to the following formula [Ref. 5] from probability 
theory 



g(v) = 4>(v x ) | J(v x /v) 



(28) 



where J (v /v) is the Jacobian and in this case equals 

X 

dv /dv. Using (25) and (27) in (28) and normalizing one 

A 

arrives at the probability that a molecular oscillator of 
mass M and resonant frequency v q in a gas at temperature T 
emits a photon at frequency v 



g( v ) 




exp 



-Me 2 

-2*Tv 0 2 




( 29 ) 



This is the Gaussian lineshape function for inhomogeneous 
broadening and has a full width at half maximum given by 
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(30) 



Av 



D 



2v 



o 




It was stated at the beginning of this section that the 
lineshape is either Lorentzian or Gaussian in the limit. 

In reality several broadening mechanisms occur simultaneous- 
ly including Stark broadening, van der Waals broadening, 
and wall collisions to name a few. In the simplified case 
where two different broadening mechanisms yield comparable 
linewidths , but which are considerably larger than all other 
broadening mechanisms, and the two mechanisms are statis- 
tically independent, the observed lineshape is equal to the 
convolution of the two individual lineshapes 



If one of the lineshapes is Gaussian due to Doppler broaden- 
ing and the other is Lorentzian due to pressure broadening, 
Eq. (31) yields a Voigt profile [Ref. 6] . The total 
linewidth in this case may be estimated according to [Ref. 7] 



Equation (31) was numerically calculated on a Hewlett- 
Packard 9810A calculator for the Lorentzian and Gaussian 
profiles in Fig. 5. The normalized result is also plotted 



g obs( v ) = ' ^')g 2 ( v ') dv 



(31) 



4v * [(4 v d ) 2 * (Av p ) 2 ] 1/2 



C32J 
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in Fig. 5. All three profiles were plotted by the 
calculator on the same scale for comparison purposes. 

The program for the Voigt profile is given in the appendix 
for the inquisitive reader. 

D. LONGITUDINAL AND TRANSVERSE MODES 

The natural resonant frequency of a laser oscillator 
cavity is given by 



cos 1 



q£m 



q + (Z + m + 1) 



/a~jt 



-)(1 



R„ 



2nL 



( 33 ) 



where q is the longitudinal mode number, Z and m are 
transverse mode numbers, R^ and R£ are the radii of curva- 
ture of the cavity mirrors, c is the speed of light, n the 
refractive index, and L is the mirror separation. It is 
easily verified from this equation that the frequency 
separation between two adjacent longitudinal cavity modes is 
given by 



v 



q + i 




(34) 



A laser which has gain over this frequency difference, i.e., 
a linewidth greater than or equal to the value of (34) , will 
support more than one longitudinal mode at a time. Thus 
a single mode device should have a short cavity length. 
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The frequency or mode separation for adjacent transverse 



modes is given by 



cos ]/ (1 - ] (1 - £— ) 



1 2 



c 



v £+m+l v &+m 



7T 



2nL 



(35) 



This quantity is considerably less than that of (34), and 
often less than the line width. The problem of higher order 
modes oscillating simultaneously is sometimes remedied by 
the use of an aperture. Bloom [Ref. 8] states that the 
lowest order mode will predominate for laser bore diameters 
less than 3.5 xoj where w is the spot size on the output 
mirror. 

The spot size for a resonator with mirrors of unequal 
radii of curvature is given by [Ref. 9] 



Dividing the numerator and denominator on the right-hand side 




r R x + R 2 - L 



(36) 



of (36) by R^“ and taking the limit as R^ 00 gives the 
spot size when the output mirror is an optical flat; 




( 37 ) 
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E. FREQUENCY PULLING 

The passive cavity linewidth for R equal to the square 
root of the product of the two mirror reflectances and a 
equal to the average distributed loss constant is given by 
[Ref. 3] 



Av 



c 



c [a - i 2/n R] 
2rrn 



c(l - R) 

2irnL 



(38) 



where the approximation holds for a = 0 and R 15 1 . The 
cavity linewidth is generally less than the Doppler broadened 
linewidth and results in hole burning in the gain curve of 
an inhomogeneously broadened system at the frequency or 
frequencies of oscillation. 

In general the passive cavity resonance v does not 
coincide with the molecular line center frequency v q . As a 
result the oscillation will occur at a frequency near v 

H 

but slightly shifted toward v Q . This is called frequency 
pulling. This frequency pulling is due to the negative 
anomalous dispersion curve which results in a shift in the 
index of refraction near the molecular resonance. The 
frequency of oscillation of the laser for a high Q cavity 
is given by 



v 




Av c 

^ V q v o^ Av 



(39) 



where Av is the total molecular linewidth due to Doppler and 
pressure broadening. 
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Since the longitudinal mode number q is immensely 
larger than the transverse mode numbers Z and m, Eq. (33) 
may be simplified to 



v 



q 




( 40 ) 



Substitution of (38) and (40) into (39) gives for the 
frequency of oscillation 



v 



= 4c_ 
2nL 



~ qc 2 (i - R) 

~a 271 
4TTn L 



v 0 c (l - 
2itnL 




(41) 



Equation (41) is an important result as it reveals the fact 
that the laser oscillation frequency may be tuned simply by 
changing the cavity length or broadening the linewidth. 

This is of fundamental importance in stabilization tech- 
niques utilizing a piezoelectric driven mirror in the cavity. 

F. INSTABILITY AND CAUSES 

Taking the derivative of (41) with respect to length 
yields 



dv 

cTT ' 




~ qc 2 (1 - R) 
- 2Trn 2 L 5 



v o c(l - R) - 
2TmL Z 



1_ 

Av 



( 42 ) 



Substituting Av c and from (58) and (40) back into (42) 
gives the change in frequency dv for a change in length dL 
as 
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( 43 ) 



dv „ ~ V q , 
dL L 





Av 

c 

Av 



The quantity in the parentheses of (43) is approximately 
v ; thus (43) is simplified to 

H 



dv _ "_^q . V q Av c 

31 L raT 



(44) 



Equation (44) shows that when Av c << Av the change in 
frequency with a change in length is large. As the cavity 
linewidth increases or the molecular linewidth decreases 
the laser oscillation frequency fluctuates less with 
corresponding length fluctuations. Normally the cavity 
linewidth is much less than the molecular linewidth; 
consequently the primary sources of instability are those 
which affect the optical length of the cavity. 

Equation (44) also yields the following approximate 
expression [Ref. 10] 



dv -dL 

v~~ L 



(45) 



Since ^ is proportional to the thermal coefficient of 
linear expansion a times the change in temperature, it 
follows from (45) that the frequency of oscillation fluctuates 
with temperature according to 
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dv 

v 



-adT 



(46) 



Equation (46) shows that an increase in temperature yields 
a decrease in frequency. Mocker [Ref. 11] gives a coefficient 
of linear expansion for low expansion quartz of 2 x 10' /°C. 
With this value of a the laser coolant must be held to 
within 1°C in order to achieve a frequency stability of one 

7 

part m 5 x 10 . 

Bridges and Patel [Ref. 12] point out that there is a 
narrowing of the Doppler linewidth with decrease in 
temperature, which is evident from (30), and consequently 
one would expect better frequency discrimination between 
competing modes. Thus if mode competition is a problem 
where an interferometric device or grating cannot be in- 
serted within the cavity, lowering the coolant temperature 

may be the solution. It is also stated in Ref. 12 that the 

-3/2 

optical gain is approximately proportional to T ' . Thus 

lower temperatures yield higher output powers as well as 
increased stability. 

Pressure fluctuations are another source of instability, 
though insignificant in sealed lasers. According to Mocker 
[Ref. 13] frequency shifts due to gas pressure changes of 
5-8 MHz per torr have been observed. This shift is toward 
the red with increasing pressure. The shift is due to a 
change in the refractive index of the gas which effectively 



33 



changes the optical length of the cavity. For a sealed 
cavity, pressure fluctuations would result from changes in 
the gas temperature especially during warm up. 

Acoustic vibration is a major source of instability. 

It also results in cavity length fluctuations. It is 
minimized by utilizing internal cavity mirrors and isolating 
the cavity as a whole from the support members with acoustic 
absorbers such as foam rubber pads. The whole laser may 
then be enclosed in an anechoic box. Significant vibration 
results from the circulating coolant water and fluctuations 
in the tap pressure. This is minimized by maintaining the 
flow rate as low as possible. Air cooled systems and motor 
driven choppers can generate copious amounts of acoustic 
noise . 

Fluctuations in the excitation current result in both 
frequency and amplitude instability. Mocker [Ref. 13] 
observed a shift of 500-900 kHz/mA shift toward the blue 
with increasing excitation. The amount of shift decreases 
with increasing excitation. The amplitude or output intensity 
is observed to have an optimum value of excitation and 
decreases with a change from this value in either direction 
[Ref. 12]. A laser which is to be amplitude and frequency 
stabilized must have a well regulated and filtered power 
supply. 

The current dependence of the oscillation frequency is 
probably due to a change in the refractive index of the active 
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medium as a result of density changes of the electrons, 
ions, and neutral molecules. The index of refraction as a 
function of the number densities of electrons (N e ) , ions (N^) , 
and molecules (N m ) is given by 



, N e 4 
_ , 1 e 

n 0 - 1 • 2 : 

me 0 (u 



2tt 7 a-N. + 2iTa N 
* i l mm 



(47) 



A reduction in the refractive index may arise from the more 
rapid creation of components such as CO and NO which have 
a smaller electric permittivity than CO^ and thus contribute 
to an increase in the frequency of oscillation. 

The total fractional change in frequency is equal to 

/ 

the sum of the fractional fluctuations [Ref. 10] , 



Av _ An AL 

v~ ~ E“ ’ 



(48) 



where An is found from (47) and AL from (45) and (46) . 

G. STABILIZATION TECHNIQUES 

As was noted earlier the output amplitude of a laser is 
primarily a function of the discharge current. Therefore, 
current regulation should eliminate a major portion of the 
amplitude fluctuations. It is also noted that due to a 
frequency-dependent gain curve the output amplitude is 
somewhat dependent on the frequency of oscillation. 
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Consequently frequency stabilization will also produce 
some amplitude stability. The following discussion of 
some of the techniques used focuses on frequency stabilization. 

In the above discussion of instabilities, it was 
observed that each fluctuation effected the optical length 
of the cavity, thereby changing the frequency. It follows 
then that frequency stability would be achieved by cavity 
length correction. This is often accomplished by mounting 
one or both morrors of the cavity on a piezoelectric crystal 
which expands or contracts in proportion to a voltage 
applied across it. An error signal is then produced by 
some method and drives the piezoelectric crystal. A 
suitable error signal must meet two basic criteria; it must 
be of a magnitude proportional to the amount of length 
correction required and it must indicate the direction of 
the required length correction. 

The various stabilization techniques may be basically 
divided into two categories, those which stabilize to a 
passive cavity resonance and those which stabilize to an 
atomic or molecular transition [Ref. 14] . 

An example of stabilizing to a passive cavity is to pass 
the laser beam through a Fabry-Perot etalon, the plate 
separation of which is modulated at an audio frequency. 

The effect of the audio frequency is to dither the cavity 
resonance frequency about the desired optical frequency, 
thus changing the transmission of the etalon. The output 
beam of the etalon is a phase and amplitude modulated signal 
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which is then converted to an electrical signal with a 
photodetector. The signal from the detector is compared 
with the audio modulating frequency in a lock-in amplifier. 
The lock-in amplifier D.C. output is the error signal. 

The disadvantage of using a passive cavity is that it is 
still subject to thermal and acoustic fluctuations. 

One example of stabilizing to a molecular resonance 
is to use the Lamb dip of the Doppler broadened lineshape. 

The Lamb dip or power dip is the result of hole burning 
in the gain curve. The laser frequency may be dithered 
about the dip and the laser output compared with the 
dithering signal to produce an error signal. This method 
has the inherent disadvantage of having a very narrow 
frequency region in which stabilization can occur. If the 
laser frequency should inadvertently drift out of the region 
of the power dip, the error signal would have the wrong 
polarity and force the laser frequency further from the 
desired resonance. 

A better method of stabilizing to a molecular resonance 
is to use an external gas-filled cell whose center absorption 
frequency is dithered about a point near the laser frequency 
by application of an oscillating magnetic or electric field. 
This type of center frequency modulation utilizes the Zeeman 
effect or Stark effect, respectively. The output of the cell 
is then phase sensitively detected as before to produce an 
error signal. 
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H. STARK CELL STABILIZATION 



If an absorbing or radiating molecule is immersed in 
a strong electric field its energy levels and consequently 
its spectrum will be changed. This is the Stark effect and 
is similar to the Zeeman effect where a magnetic field is 
used as the perturbation rather than an electric field. 

The application of the electric field removes the degeneracy 
of the vibrational -rotational levels and results in a splitting 
of the levels into several other levels. Kelly, Francke, and 
Feld [Ref. 15] give the relevant energy level diagram for 
absorption of several CO^ laser lines by NF^D and show the 
splitting of the vibrational -rotational levels. 

There are many gases whose absorption peaks can be 
brought into coincidence with CC^ laser lines by the Stark 
effect. Martin, Corcoran, and Smith [Ref. 16] give a table 
of different gases and the laser lines absorbed. The 
remainder of the discussion, however, will focus on the 
absorption of the P(20) C0 9 laser line in the 10.4 urn band 
by NH 2 D. Brewer, Kelly, and Javan [Ref. 17] found that this 
coincidence could be obtained at field strengths of 3298 
V/cm, 4404 V/cm, and 6604 V/cm. 

The D.C. Stark effect can be represented from perturba- 
tion theory to first order by 

E m = E m - ’ MV (49) 
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where E^ is the unperturbed energy of state m, is the 
unperturbed wave function, p the molecular dipole moment, 
and e s the Stark field operator [Ref. 19] . 

The transition frequency between two Stark shifted 
levels is 



v 

mn 




(50) 



Substituting a time -dependent Stark field of the form 
£ s (t) = e s cos cot into (49) yields 

E m (t) = E m ' < f' m l P ’ c slV cos Mt < 51 ) 



where the cosine has been brought outside the integral 
since the integration is over distance. It is also assumed 
in (51) that to << v mn ; otherwise the energy levels could 
not follow the field. Equation (50) then becomes 
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Equation (52) can be rewritten as 

v (t) = v + Av ’ cos cot . (53) 

mn o 
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Equation (53) is an important result which will be used in 
a later analysis. It shows that the application of an 
oscillating field results in a frequency modulation about 
a center frequency v . Claspy and Pao [Ref. 18] solve the 
equation of motion using the density matrix approach and 
show that for Stark field frequencies less than the homo- 
geneous linewidth Av , the absorption coefficient is a 
single Lorentzian line with a time varying center frequency 
moving synchronously with the Stark field. The result is 
similar to (53) . 

Nussmeier and Abrams [Ref. 19] give the frequency dif- 
ference between the P(20) 10.6 um line center and the 
absorption center of NH 2 D as a function of applied field as 

6 = 2042 + 0 . 143 | M | E (54) 



where 6 is in MHz, M is the z component of the rotational 
quantum number J and can take on integer values from -4 to 
4, and E is the applied electric field in V/cm. Equation 
C54) shows that for M = ±4 an electric field of 3570 V/cm 
is required to bring the P(20) line into resonance with 
NH 2 D. Equation (54) differs only slightly from the results 
found by Brewer, Kelly, and Javan [Ref. 17]. The quantity 
0 . 143 | M | in (54) is called the Stark tuning rate and is 
reported only slightly differently in the literature. For 
example, Tangonan and Abrams [Ref. 20] report a tuning 
rate of NH 2 D and the P(20) line of .564 MHz cm/V for M = 4 
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which gives 0.141 |M|, Johnston and Melville [Ref. 21] 
report a tuning rate of .568 MHz cm/V. 

A proposed stabilization scheme is shown in Fig. 6, 
which is a slightly modified version of that given in 
Ref. 19. The Stark cell D.C. power supply biases one plate 
of the cell to bring the P(20) line into resonance with 
the NH^D gas. This occurs at a field strength of 3570 V/cm. 

The other plate is driven with an AC signal producing a 
field strength of approximately 150 V/ C m at about 10 KHz 
[Ref. 19]. The cell is filled with a 1:1 mixture of ND^ and 
NH^ at approximately 2 torr. The exchange occurs rapidly 
to produce N^D [Ref. 21]. 

A beam splitter reduces the radiant flux incident on 
the Stark cell below lW/cm" to avoid saturation effects. 

The transmitted intensity is then detected with an SBN pyro- 
electric detector or HgCdTe detector. The detector output 
is fed to the "signal in" on the lock-in amplifier. The 
audio oscillator signal is fed to the "reference in" on 
the lock-in amplifier. The lock-in determines the phase 
of the detector signal with respect to the audio oscillator 
input. The polarity of the D.C. lock-in output signal is 
determined by the phase of the detector signal. The ampli- 
tude of the D.C. output is proportional to the amplitude 
of the detector signal. The D.C. lock-in output signal 
then is the error signal. 

The error signal is integrated over time in the integrator. 
The integrator output drives the programmable D.C. power 
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FIGURE 



supply which in turn drives the piezoelectric-mounted laser 
cavity mirror. Thus the loop is completed. The exact 
details of error signal production and the integrator are 
covered in the experimental part as a result of a computer 
analysis . 
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III. EXPERIMENTAL 



The following discussion of experimental procedure and 
results is specific to the GTE Sylvania Model 941 C0 2 Laser 
whose characteristics are given in Appendix A and to the 
Stark cell whose characteristics are given in Appendix B. 
Unless otherwise indicated, wherever the phrases "the laser" 
or "the Stark cell" appear, it is understood that they 
refer to the laser and Stark cell described in Appendices 
A and B . 

A. LINESHAPES 

An attempt to determine the lineshapes of the laser 
transitions experimentally was made by plotting relative 
intensity and mode as a function of piezoelectric voltage. 
The experimental arrangement is schematically represented 
in Fig. 7. The x-axis of the x-y plotter sampled the 
piezoelectric voltage from terminals 8 and 10 of the PZ 
power supply through a voltage divider. The PZ power supply 
is a Kepco Model ABC regulated D.C. supply with a range of 
0 to 1500 volts. It was operated in the voltage program- 
mable mode with a positive ground. The programming voltage 
was a regulated 0 to 25 volt D.C. supply. The battery box 
and reversing switch were added to obtain the proper biasing 
for tuning over the full range of the Kepco with a minimum 
of switching. 
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FIGURE 



The y-axis of the x-y plotter sampled the relative 
intensity of the chopped laser output via a pyroelectric 
detector, preamplifier, and a filtered rectifier. The 
longitudinal mode or vibrational - rotational transition was 
continuously monitored using an Optical Engineering C0 2 
Spectrum Analyzer. Frequency could not be monitored due 
to the lack of precision, but the longitudinal mode could 
be identified. Typical results are given in Fig. 8 for an 
increasing PZ voltage and in Fig. 9 for a decreasing PZ 
voltage . 

It is noted that the plots are not identical but a 
hysteresis effect exists [Ref. 22]. The P(20) line is 
shifted in the direction of the voltage scan. Various scan 
rates were tried in order to ascertain that the effect was 
not due to the response of the instruments. Nine to ten 
different lines were observed over the tuning range. It 
is also noted from Figs. 8 and 9 that the series of lines 
begins to repeat itself after approximately 1.4 kV. . Since 
the resonance condition requires the cavity length be an 
integer number of half wavelengths, this voltage change 
corresponds to a length change of half a wavelength and the 
piezoelectric ^ parameter may be estimated as 

AL s 1 / 2 ( 10.6 pm) = 3 g um / kV (55) 

Av 1 . 4 kV 



46 



INTENSITY AND MODE AS FUNCTION OF 
INCREASING PZ VOLTAGE 
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PZ VOLTAGE 



It is noted that this value corresponds with the value 
given by Sylvania. This calculation assumes that the 
piezoelectric crystal has a linear response to voltage. 

From (55) and Figs. 8 and 9 the apparent linewidth 

of the P(20) transition is estimated to be 55 MHz which 

is close to the Doppler linewidth given by Sylvania. The 

cavity linewidth calculated from (38) with a = 0 is 19 MHz. 

The pressure broadened linewidth is 100 MHz and the 

Doppler linewidth is 50 MHz. The total linewidth using 

(32) is 112 MHz. Taking the differential of (39) and using 

19 MHz for Av c , 112 MHz for Av , and 55 MHz for the change 

in passive cavity resonance dv , gives a change in oscilla- 

4 

tion frequency dv of 46 MHz, which would be closer to 
the actual linewidth if the above numbers are correct. 

It is concluded that the PZ voltage may be used to 
estimate the laser oscillation frequency; however, the 
actual frequency change is only slightly less than the change 
indicated by the PZ voltage due to the frequency pulling 
effect. A more precise means of measuring the actual 
frequency of oscillation would be required to verify (39) 
and determine the average distributed loss constant a. 

For gross changes in the frequency of oscillation, using the 
PZ voltage alone will result in insignificant error. To 
determine small changes such as linewidths a more precise 
method of measuring frequency is required. It would also be 
necessary to insert a prism or grating in the cavity to 
prohibit the neighboring modes from oscillating until the 
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half power point was reached. From Figs. 8 and 9 it is 
observed that the half power points are never reached on the 
P C 2 0 ) line before the neighboring line begins to oscillate. 
There is considerable overlapping of the lines. Mode compe- 
tition is also quite prevalent in a CO 2 laser. Thus it is 
impossible to measure linewidths in the manner attempted above. 

B. MODES OF OPERATION 

The laser's cavity length is .45 m. This gives from (34) 
a mode separation of 333 MHz which is considerably larger 
than the total linewidth of 112 MHz. Consequently the laser 
should normally oscillate in a single longitudinal mode, 
which in fact it did. Occasionally two lines were observed 
to oscillate simultaneously and on one occasion four lines 
oscillated simultaneously at a coolant temperature of about 
18°C; the P ( 2 2 ) , R(16) R(24) , and R(30) lines. Also a very 
faint P (14) line of the 9.4 ym band was observed to oscillate 
simultaneously with P(22) on the 10.4ym band at 19-20°C 
coolant temperature out. 

The laser utilizes a dielectric flat for the output 
mirror and a 3 m totally reflecting mirror mounted on a 
piezoelectric crystal. Using (35) gives a transverse mode 
separation of 42 MHz, which is less than the total line- 
width. From (37) the spot size on the output mirror is 
1.9 mm. The laser's bore diameter is 6 mm. Thus from the 
discussion preceding Eq. (36) the laser is only marginally 
diffraction limited to the lowest order transverse mode. On 
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one occasion the TEM^q/TEMq^, the "doughnut" mode, was 
observed. On another occasion efforts to reproduce this 
mode yielded only the TEMqq mode throughout the full 
range of the PZ voltages. 

C. INSTABILITIES 

Experimentation began with minimization of 60 Hz 
noise. The voltages on the ground terminals of all 
115 VAC receptacles in the lab were measured with respect 
to a drainpipe ground. One receptacle ground was found to 
have a 52 mV 120 Hz signal. Its ground strap was isolated 
from the grounds of the other receptacles in order to 
minimize the signal on them. The receptacle with the 
lowest ground signal (5 mV peak to peak at 180 Hz) was 
used as the source for all measuring instruments and power 
supplies. Noisy receptacles were used for vacuum equipment, 
chopper, and heat lamps. Ground loop signals were eliminated 
only by trial and error. 

The laser power supply, the GTE Sylvania model 941P, 
is unregulated and minimally filtered. It was found to have 
a one to ten per cent current ripple depending on the dis- 
charge current. This was entirely unsuitable for stabiliza- 
tion of the laser as the current ripple showed up in an 
amplitude modulated laser output. An L-section filter 
consisting of a . 5 pf capacitor and a 45 H choke (which was 
actually the secondary of a 13 kV transformer) was assembled 
and connected to the high voltage output at the laser power 
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supply. Details of the filter are given in Appendix E. The 
filter reduced the current ripple and consequently the laser 
output fluctuation by a factor of 50 to approximately .2%. 
Since a 10% Stark cell modulation was expected, this filter 
was considered to be suitable for short-term stability 
purposes. For long-term stability, a well regulated power 
supply was still required. Plots of laser intensity as a 
function of excitation current were made using an x-y 
plotter. The x-axis input was the voltage drop produced 
by the excitation current across a small resistance element. 
Typical results are given in Fig. 10. The relative intensity 
variation is approximately 14%. 

The experimental setup for determining the effect of 
cooling water temperature on the intensity and longitudinal 
mode is schematically represented in Fig. 11. The coolant 
output temperature was sampled using a copper-cons tantan 
thermocouple referenced to an ice bath. The voltage drop 
across the thermocouple was measured with a nanovoltmeter 
on the 3 millivolt scale. The output of the nanovoltmeter 
drove the x-axis of an x-y plotter. The y-axis sampled 
the chopped laser intensity via a pyroelectric detector, 
preamplifier and filtered rectifier. The coolant tempera- 
ture was slowly increased over a period of about one hour 
from 7° to 40°C by circulating the coolant through an ice 
bath and warming the ice bath with a heating element. The 
longitudinal mode was continuously monitored using a spectrum 
analyzer. Mode changes were noted and recorded on the plot 
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INTENSITY & MODE AS A FUNCTION OF 
TUBE CURRENT 
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TUBE CURRENT 



INTENSITY VS TEMPERATURE EXPERIMENT 
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ICE BATH FIGURE 11 



after each run. A table of temperature as a function of 
thermocouple voltage is given in Appendix F. 

Figure 12 is a plot with the PZ voltage held constant. 
Note that the series of lines from P (18) to R(22) is repeated 
in identical order; however, there is little similarity to 
the intensity profiles. The P(18) line is repeated twice 
over a temperature range of 8.8 to 37.9°C. Since the 
resonance condition requires the cavity length to be 
exactly an integer number of half wavelengths of the oscilla- 
tion frequency, this temperature range corresponds to a 
thermal expansion of two half wavelengths or one wavelength 
of the P(18) transition. Thus the thermal coefficient may 
be determined from Eqs . (45) and (46) as 



1 dL = 1 (10.57 x 10‘ 6 ) 

" L 3T 7TT 29.'1 U C 



8 x 10' 7 /C° . 



(56) 



Note that this calculation assumes that a temperature change 
in coolant out corresponds to an equal temperature change 
in the laser tube wall. See Appendix G for proof of this. 

The plot shows an overall decrease in intensity output 
with increasing temperature, by a factor of about two in a 
temperature range of 33C°. It is also obvious that the P(20) 
transition is the strongest in the group. P(22), R(16) and 
R(24) form a competing trio. Mode competition is almost 
always present between these lines. All three or pairs of 
these lines are often observed oscillating simultaneously. 
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INTENSITY & MODE AS FUNCTION OF 
INCREASING TEMPERATURE 

(COOLING WATER OUT) 
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thermocouple voltage 



It would therefore be wise to avoid stabilizing on any one 
of these three lines. 

Figure 13 shows the maximum intensity of the P(20) line 

as a function of temperature. The intensity decreases 

almost linearly with increasing temperature. The region 

between 7 to 27°c shows mode competition. It is also noted 

from Fig. 13 that as temperature increases or the cavity 

expands, the PZ voltage had to be increased in order to 

maintain the P(20) line. It follows that an increasing PZ 

voltage corresponds to a decreasing cavity length and the 
, AL 

sign on the parameter given by (55) should be negative. 

An experiment was performed to determine the effect 
of coolant flow rate fluctuations on the laser cavity. 

The flow rate was adjusted within the limits set by Sylvania 
from 11 ml/s to 30 ml/s at a constant temperature input. 

The temperature output changed as expected from 23.8 to 
22.8°C. Using (56) this one-degree temperature change 
corresponds to a length change of .36 pm. In this experiment 
the PZ voltage was adjusted to maintain maximum intensity 
on the P ( 2 0 ) line, and was adjusted from .42 kV to .24 kV. 
Using (55) this voltage change corresponds to a length change 
of .68 pm which is almost twice the length change accountable 
for temperature. It follows that a flow rate fluctuation of 
19 ml/s results in a length fluctuation of approximately 
.32 pm. From (45) this corresponds to a frequency fluctua- 
tion of 22.6 MHz. 
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INTENSITY & PZ VOLTAGE OF P(20) MODE 

AS FUNCTION OF INCREASING TEMPERATURE OF 
COOLING WATER OUT 
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THERMOCOUPLE VOLTAGE 



It is concluded from these experiments that temperature 
and flow rate fluctuations are the dominant causes of 
instability in the laser. Current fluctuations can be 
easily remedied with a filtered current regulated power 
supply. If a circulation pump is used for cooling, the 
system should be well designed to minimize pressure and 
flow rate fluctuations . Experiments could not be carried 
out to determine the effects of acoustic vibration due to a 
lack of a method for precisely measuring small, rapid 
fluctuations in frequency. 

D. ERROR SIGNAL ANALYSIS FOR A PYROELECTRIC DETECTOR 

Drawing an analogy with electrical communication theory 
it is observed that the stabilization technique is similar 
to frequency demodulation. The input laser beam is 
frequency modulated with noise. It is desired to have some 
sort of frequency discriminator which will sense the shift 
in frequency and produce an output voltage linearly pro- 
portional to the shift. This is precisely what FM demodu- 
lation accomplishes [Ref. 23] • The stabilization servo 
loop is completed when the demodulation information is fed 
back to the laser to make corrections to minimize the noise 
modulation . 

The exact details of producing an error signal were 
not clear in the literature. Questions which needed to be 
resolved were: is the error signal produced by slope or 

peak detection, how does dithering produce an error signal, 
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is the composite detector signal or one of its harmonics 
monitored, and which harmonic should be monitored? 

In order to search for a suitable error signal we 
begin by assuming the transmitted intensity of the cell 
of length L is of the form 

I(v) = I Q (v) exp [ - y (v) L] (57) 

where I Q C v ) is the input intensity, y(v) is the absorption 
coefficient, and I(v) is the output intensity as a function 
of frequency. It is assumed that I (v) is a narrow spectral 
laser line of half width much less than the cell absorption 
linewidth. Thus I (v) may be approximated by a Dirac delta 
function, <5(x). If the frequency of the laser is allowed to 
vary randomly with time over the linewidth of the absorption 
cell, and the linewidth of the absorption cell is assumed to 
be much less than the laser's molecular linewidth so that the 
laser intensity does not change appreciably over the cell 
linewidth, then (57) becomes 

I v (t) = n 0 <5(x(t) - v)e' Y(v)L dv = I o exp [-y(x(t))L] (58) 

where x(t) is a random variable of time and takes on values 
of frequency within the cell's absorption linewidth. The 
subscript v on I denotes that I is approximated by a delta 
function, i.e., all the intensity is concentrated at a 
single frequency. 
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Equation (58) shows that the frequency modulation x(t) 
has been converted to an intensity modulation I(t) by simply 
passing the laser beam through a selective absorber or 
optical notch filter. It is possible to derive an error 
signal from I(t) by stabilizing about a point on the slope 
of y(v) . This is called slope detection, analogous to slope 
detection in radio FM demodulation. Two problems arise with 
slope detection; spurious intensity variations in the laser 
beam also show up in the error signal and the range of the 
linear slope is quite small [Ref. 23] . To partially over- 
come these problems, it is desirable to obtain an error 
signal which is similar to the frequency to voltage charac- 
teristic of a common FM balanced discriminator, where the 
voltage goes to zero at the center frequency, increases 
positively on one side of the center frequency, and increases 
negatively on the other side. 

Consider the effect of modulating the center frequency 
of the absorption cell with an audio signal of frequency to 
according to (53) , and using phase-sensitive detection in 
order to determine zero time. This is schematically repre- 
sented in Fig. 14. The modulating signal will dither the 
absorption profile about v Q at frequency to. The amplitude 
of the dither is Av ’ . It is assumed that the dither is 
much more rapid than the random variation of the laser 
frequency v. As will become more apparent later the reason 
for the assumption is to detect x(t) . Only fluctuations 
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FIGURE 14 
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slower than w can be detected. The end result of the 
analysis to follow is that x(t) amplitude modulates a 
"carrier frequency" w. It would be desirable to make co 
extremely large in order to detect as much noise as possible. 
However, as Claspy and Pao [Ref. 18] show, the upper limit 
for a) is the homogeneous linewidth. This is reasonable if 
one considers that from (18) the linewidth was derived from 
the decay time. It is now observed that the decay time is 
also the absorption time. Thus if a molecule does not 
remain in resonance with the passing radiation for the 
duration of the absorption time, the molecule will not 
absorb. In other words, the Stark cell is bandwidth 
limited by the homogeneous linewidth. 

The detector used in this part of the analysis has an 
output voltage proportional to the time rate of change of 
the incident intensity. An example of such a detector is a 
pyroelectric. This time-differentiating detector will be 
referred to as "detector A" in this analysis. 

The absorption coefficient of the cell is assumed to 
have a Lorentzian shape from (19) . 



where S is the line strength and is a constant. The absorp- 
tion coefficient is given a time dependence by oscillating 
the center resonant frequency according to (53). 



y(v) = Sg(v) 



SAv 



(59) 




(v 
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Y(v,t) 



SAv 



(60) 



2tt [(y~) ^ + (v - v q - Av 1 cos »t) 2 ] 

The difference in frequency between the laser beam and the 
average center resonant frequency v Q is defined as 
5 = v - v Q , where the symbol used is not to be confused 
with the Dirac delta used earlier. Equation (60) then 
becomes 



Y(v,t) 



SAv 

2 it [(^-) 2 + (5 - Av' cos ut) 2 ] 



(61) 



Substituting (61) into (58) the intensity out as a function 
of time is given by 



I(t) = I Q exp 



SAvL 

2tt [(^O 2 + (6 - Av 



cos 




(62) 



Detector A output as a function of time is a voltage 
proportional to 

f(t) = = -I q SL exp {-SL g(v,t)} (63) 

As was assumed earlier the oscillation frequency is much 
greater than the random variation of the laser frequency v. 



64 



Thus 6 is approximately constant and the derivative with 
respect to time of (62) may be taken holding it as such. 
Taking the derivative of (59) , expanding all squared terms 
and using common trigonometric identities yields 



g dg(v,t) 

cTt 



-SAvAv ' oo(6 sin tot - sin 2wt) 
tTLA+b-cos tot + C • cos 2 cot + U-cos boot + E'*cos 4cotj 



where 



A = 6 4 + 36 2 (Av') 2 + -g-(Av') 4 + (j^-) 2 26 2 + (^-) 2 (Av') 2 + ( 
B = - 46 °Av ' - 36(Av’)° - 46Av'(y^-) 2 

C = 36 2 (Av ' ) 2 + j(Av') 4 + (Av*) 2 (^-) 2 (fi4) 

D = -6 (Av ' ) ^ 

E = j-CAv') 4 



Equation (63) can be simplified if y(v,t)L << 1. The 
exponential then is approximately one and f(t) is pro- 
portional to (64) . Note that S could be determined 
experimentally by applying a D.C. field and measuring the 
transmitted intensity when v = v Q and when |v - v | >> Av. 

Equation (64) was plotted on a HP 9810A calculator 
for various values of 6 and Av ' . The full width at half 
maximum Av and the period of the fundamental oscillation T 
were arbitrarily set to one in the program; thus go the 



Av-, 4 

T~ J 
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fundamental frequency is 2tt. 5 is varied between -1.0 and 
1.0 and Av ’ between 0.1 and 1.0. The results are plotted 
in Figs. 15, 16, 17, and 18. Figure 15 is for 5=0 and 
Av ’ = .1 and .7. Note the frequency is twice the funda- 
mental. This would be expected for a laser line oscillating 
symmetrically about the peak of a transfer characteristic. 

From (64) it is also noted that when 5 = 0 no fundamental 
or first harmonic exists. Figure 15 also reveals the non- 
linearity produced by overdriving with a large Av ' . Though 
the output signal is greater, higher order harmonics are 
obviously present, producing a distorted sin 2wt. 

Figure 16 is for 5 = .2 and Av f = .1 and .7. It is 
obvious that some first harmonic is present. It is noted 
that when Av ' > 5 the dither carries the line to both sides 
of the y(v) peak, thus producing a large amount of second 
harmonic but because the dither is unsymmetrical about 
the peak of y(v) a small amount of first harmonic is 
generated. When 5 > Av' second harmonic is virtually 
nonexistent . 

Figure 17 is for 5 = .5 and Av ’ = .1 and .7. This is 
the region of slope detection. A considerable amount of first 
harmonic is obviously present. Very little second harmonic 
is present even when Av ' is greater than 5. Figure 18 shows 
a 180° phase shift in going from a positive 5 to a negative 5. 

It is concluded from these four plots that the basic 
requirements for generating an error signal are present; the 
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DETECTOR A OUTPUT as a function of time 

FOR TWO PERIODS OF FUNDAMENTAL, §=0 
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FIGURE 15 
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FIGURE 16 
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FIGURE 17 



DETECTOR A OUTPUT 

AS A FUNCTION OF TIME 
FOR TWO PERIODS OF 
FUNDAMENTAL , 




70 



FIGURE 18 



amplitude of the signal indicates the magnitude of 6 and the 
phase of the signal indicates the sign of 6. Thus the 



distance of the laser frequency from the reference v Q is 

known and the direction from v is known. It is still 

o 

desired, however, to obtain the typical FM frequency dis- 
criminant and to determine which harmonic to monitor or 
if the composite signal should be monitored. 

It is revealing to decompose the signal of (64) into 
all the harmonics which make up the signal and examine 
each harmonic singly. Any periodic signal of fundamental 
period T which is bounded and real may be expanded into a 
Fourier series according to 



f(t) 




+ y (a 

n 

n=l 



cos 



n2irt j. 
T 



b 



n 



sin 



n2irt~, 
T J 



where 



T 

a n = % f f(t) cos dt , n = 0, 1, 2 ... 

1 o 

b - £ / f(t) sin dt , n = 1, 2, ... 

For the complex signal given by (64) the integration 
cannot be performed analytically; however, it can be per- 
formed numerically on a computer. It is known that for an 
odd function of t, i.e., f(-t) = -f(t), which is what the 
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plots show (64) to be, all the a^'s will integrate to zero 
over one period. Thus only the b^'s need be calculated. 
For a fundamental period of T = 1, Eq. (65) reduces to 
solving the problem 

00 

f(t) = l b sin (n2irt) 
n=l 



1 

b n = 2 / f(t) sin (n2TTt) dt n = 1, 2, ... (66) 

o 

for f(t) proportional to (64). 

A program for the HP 9810 was written to do the numerical 

integration of (66) up to n = 10 for various values of 6 

and Av ' and print the results. It was concluded that the 

Fourier coefficients b n for harmonics greater than n = 3 

were too small to be of any value in stabilization. A 

program was then written to calculate the Fourier coefficient 

b n for a particular n preset into the calculator for 200 

different values of 6 from -1.0 to 1.0 and plot the results. 

As before Av and T were arbitrarily set to one. Av 1 and 

o 

5 have the same units as Av and are to be viewed as a 
fraction of the linewidth Av. Av may be scaled up to any 
desired linewidth; 5 and Av ' are then scaled up with the same 
scaling factor. This program is given in Appendix D. 

Plots of the first three Fourier b coefficients of 
(64) as a function of <$ , with S = 1 , are shown in Figs. 19, 
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20, 21, and 22. Figures 19 and 20 are the typical shape 
of an FM discriminant. These plots represent the relative 
magnitude of the D.C. error signal produced by a lock-in 
amplifier when only a single harmonic is phase sensitively 
detected with the oscillator signal to the cell as the 
reference. Other harmonics could be filtered out using a 
narrow band pass filter. Figures 19 and 20 show that the 
first harmonic should be monitored. Note the linear region 
near 6=0. When the laser frequency is at v q the error 
signal is zero. As the laser frequency drifts to the right 
of v q for a positive 6 the error signal increases linear 
and negative. A drift to the left of v q would give a posi- 
tive error signal. This is called peak detection as the 
laser frequency is stabilized about the peak of the absorp- 
tion profile. 

Figure 20 shows that a large dither amplitude produces 
a suitable error signal. Note that the peaks are further 
away from center than for a small dither amplitude. Figure 
21 shows that the second harmonic is unsuitable for pro- 
ducing an error signal. However, it might be used in 
conjunction with the first harmonic for making a gross 
correction to the laser cavity. Note that if the laser 
frequency drifts beyond the peaks of the first harmonic 
signal, the second harmonic goes negative indicating that 
the range of stability has been exceeded and a gross 
correction to the laser cavity in the direction indicated 
by the sign of the first harmonic is required in order to 
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FIGURE 19 
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FIGURE 21 
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FIGURE 22 



bring the laser frequency back into the linear portion of 
the discriminant. 

Figure 22 shows that the third harmonic may be used 
to produce an error signal if the dither amplitude is 
sufficiently large. 

It is concluded that an error signal may be produced 
by phase sensitive detection of the fundamental or first 
harmonic of the detector signal. The lock-in amplifier 
establishes zero time from the reference signal. The 
phase of the detector harmonic will either be 0 or 180°. 

The polarity of the error signal is determined by the 
phase of the harmonic being monitored. The amplitude of 
the error signal is proportional to the amplitude of the 
harmonic, and for the first and third harmonics the error 
signal is linearly proportional to the frequency difference 
6 in the region near 6=0. 

E. ERROR SIGNAL ANALYSIS FOR A HgCdTe DETECTOR 

Consider now using a nondifferentiating detector. The 
detector has an output voltage proportional to the absolute 
incident intensity. Such a detector would be a HgCdTe 
detector. This type of detector will be referred to as 
"detector B." With detector B in the loop the detector 
signal will be proportional to I(t) as given in (58). 
Equation (58) can be approximated to first order according 
to e X = 1 + x for x small. Thus the detector output voltage 
is proportional to (61) and should contain a large D.C. 
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term due to expanding the exponential. A plot of (61) as 

a function of time showed that this was indeed the case. 

The function is also an even function, i.e., f(-t) = f (t) ; 

thus expansion is a Fourier cosine series and only the a 's 

n 

in (65) need be calculated. 

Plotting the coefficients of the harmonics as a function 
of 6 as before should show whether or not detector B 
yields a suitable error signal. Plots for the first, 
second, and third harmonics are given in Figs. 23, 24, and 
25. The results are similar to detector A. Phase sensitive 
detection of the first or third harmonic would yield an 
error signal. The second harmonic is unsuitable for error 
signal production. 

F. THE SERVO LOOP INTEGRATOR 

It is desired to maintain the P(20) line in coincidence 
with the center frequency of a Stark shifted NH ? D transition 
as accurately as possible. When the two transitions are in 
coincidence the error signal is zero. Using Fig. 19 as 
an example, consider a slight temperature increase in the 
cooling water. The laser cavity will expand and pull the 
oscillation to a lower frequency according to (45) . This 
results in a negative o and a positive error signal pro- 
portional to <5 (assuming <5 > -.2 in Fig. 19). Feeding the 
error signal directly to the PZ supply then causes the 
supply to increase the PZ voltage and make a cavity length 
contraction correction. But as the cavity length returns 
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FIGURE 23 
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FIGURE 24 
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FIGURE 25 



to the original length the error signal decreases until an 
equilibrium point is reached at some finite nonzero 6, such 
that the frequency difference 6 produced by the PZ voltage 
produces just the right error signal to produce that PZ 
voltage. It is obvious that 6 could never return to zero 
as the error signal would be zero indicating no length 
correction . 

A method is needed to maintain the required PZ voltage 
as the error signal returns to zero. This could be accom- 
plished by integrating in time the error signal produced 
by the lock-in and producing a new error signal to drive 
the PZ supply which is proportional to the total time 
integral of the first error signal. To the author's know- 
ledge, the use of an integrator in the servo loop has not 
been mentioned before in the literature. 

The integrator is placed between the lock-in amplifier 
and the PZ power supply as indicated in Fig. 6. Now as 
the error signal begins to return toward zero, the integrator 
continues to increase its output until the error signal is 
actually zero. At o = 0 the integrator maintains a 
constant D.C. output indefinitely or until a nonzero error 
signal indicates a need to change. If the error signal 
again went positive the integrator output increases still 
further, until a return to zero of the error signal indi- 
cates for the integrator to maintain a constant output at 
the new voltage reached. 
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If the error signal went negative the integrator would 
lower its output and continue to do so until the error 
signal actually returned to zero. 

The ideal integrator would respond to both rapid and 
slow changes of the error signal. It would integrate 
rapidly enough so as to stabilize against acoustic vibrations 
yet also be able to accommodate slow changes in tempera- 
ture. It should be designed to stabilize over the maximum 
temperature change expected and maximum pressure change 
expected. Such a device would produce a highly stable 
system suitable for operation in the field or in the 
laboratory . 

G. THE STARK CELL 

Stabilization was never achieved due to numerous design 
problems in the Stark cell. Six different cells were made. 
These are discussed below. 

1 . Stark Cell One 

The first Stark cell had stainless steel plates 
10 cm long and 1.5 cm wide glued with plastic spacers into 
a plastic tube 1.7 cm I.D. Stainless leads fastened to 
the plates were sealed with glue through holes in the side 
of the tube. The cell was made to accommodate one -inch 
NaCl windows as were all the cells. The plate material was 
too thin and as the glue hardened on the leads, the plates 
warped. 
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2. Stark Cell Two 

This was made with heavy stainless plate 10 cm 
long and 3 cm wide with a plate separation of 1.4 mm. 

The maximum observed D.C. modulation was 20 % at a D.C. 
voltage of 490 volts and a pressure of 4.7 torr. The 
electric field was perpendicular to the laser polarization. 
The problem was that breakdown occurred at 500 volts. 

Also diffraction was a problem with the narrow plate 
separation. 

3 . Stark Cell Three 

This was made with plates 20 cm long, 3.5 cm wide 
and about 3 mm thick. Care was taken to round the sides 
and corners of the plates. The plates were separated to 
4 mm in an effort to eliminate breakdown and minimize 
diffraction. The cell had a gross air leak which could 
only be remedied by redesign. The cause of the leak was 
the rubber "0" rings. 

4 . Stark Cell Four 

This was the same as above with slightly modified 
end window mounts. Slight modulation was obtained at 1400 
volts at a pressure of 47 torr with the E field perpendicular 
to P. The intensity through the cell was .12 W. Breakdown 
also occurred at this potential. At lower pressures 
breakdown occurred at lower potentials . 

It was then pointed out that in a low pressure gas 
electricity does not always take the shortest path. In 
fact, it was observed in this cell that at certain pressures 
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the discharge occurred between the outside two surfaces of 
the plates where the distance was greater rather than 
between the two inside surfaces. It is a fact that break- 
down potential is not a function of distance but oi 
pressure times distance. A plot of breakdown potential as 
a function of p x d is a skewed parabola with a minimum 
at some value of p x d. This value is called the pd 
minimum or Paschen minimum. To the right of the pd minimum 
breakdown potential increases with increasing pressure and 
distance as would normally be expected. But to the left 
of the pd minimum the breakdown potential decreases with 
increasing pressure and distance. Consequently if one is 
unfortunate in his choice of pressure and distances such 
that he is operating to the left of the pd minimum, increasing 
the plate separation will only make matters worse. 

5 . Stark Cell Five 

The above cell was then disassembled and mylar was 
glued to the inside surfaces of the plates. All other 
exposed metal surfaces including the lead wires were 
coated with RTV silicone rubber. As the cell was placed 
under vacuum the glue and RTV began to outgas forming gas 
pockets on the plates and in the RTV. The required 1500 
volt potential placed on the plates punched through the 
RTV and still resulted in breakdown. 

The problem was reevaluated. The operating pressure 
could not be changed enough to solve the breakdown problem 
as the Stark effect ceases to give good modulation at too 
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high or too low a pressure [Ref. 18]. Two physical 
characteristics were desired which appeared to be mutually 
exclusive: a wide channel between the plates for ease of 

alignment and a minimum of diffraction, and the elimination 
of the breakdown problem. 

6 . Stark Cell Six 

A completely new design approach was taken: to 

effectively place all metal surfaces outside of the cell, 
isolated from the low pressure gas filled channel. A 
3 mm by 10 mm wide channel was cut in a solid 22 cm long 
block of plastic. The plates were embedded in the plastic 
on either side of the 3 mm width. The design called for 
1 mm of plastic between each plate and the channel, resulting 
in a plate separation of 5 mm. Due to crating of the 1 mm 
thick wall in the machining process, a film of mylar was 
added to the underside of each plate when they were embedded 
in the plastic. This resulted in an uneven plate separation 
of approximately 7 mm. The cell was tested to 3000 volts 
with no breakdown. It also held a good vacuum. However, 
it is anticipated that the uneven plate separation may 
result in a broadened or distorted absorption profile. 

Testing was not carried out further due to a permanent 
failure of the laser and its power supply. 

This design appears to be the best approach. With this 
plate separation of 7 mm, wider plates should be used to 
minimize fringe effects of the electric field. Also the 
fabrication process needs to be modified so that the plates 
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are parallel, but sealed in plastic with no crazing in the 
high field areas. The plates should be thick so that the 
corners and edges can be rounded to minimize possible 
breakdown of the plastic. One solution might be to first 
cement a heavy plate to a block of plastic then machine the 
channel and cement the two halves together. 
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IV. CONCLUSION 



The instabilities of the Sylvania C0 2 Model 941 laser 
have been studied and the laser determined to be suitable 
for Stark cell stabilization. It was found that the exist- 
ing unregulated, unfiltered laser power supply was unsuitable 
for a stabilized system and a new well regulated power 
supply needs to be purchased. Laser characteristics 
required to build a stabilized system such as piezoelectric 
parameter and thermal expansion coefficient have been 
measured. The optimum discharge current was found to be 
7.5 mA. All the physical and operating characteristics of 
the laser pertinent to the design of a stabilized system 
have been studied and listed. 

A computer analysis has been carried out to demonstrate 
the feasibility and production of an error signal with 
existing equipment. Both a Pyroelectric and HgCdTe detector 
are suitable for error signal production by phase sensitive 
detection of the first harmonic. A secondary result of the 
computer analysis was the demonstration of the need for a 
servo loop integrator. The purpose and operation of such 
an integrator was discussed. Theoretically the integrator 
should provide increased stability over that indicated 
thus far in the literature. 

A Stark cell with no known major design flaws has been 
constructed and tested to 3000 volts with no breakdown. It 
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torr. It is expected to 



is vacuum tight to at least 10 
be successfully used in the stabilization system. 

All preliminary engineering studies have been carried 
out for assembling a Stark stabilized CO^ laser system 
operating on the P(20) mode. This system will be opera- 
tional in both the laboratory and the field. 
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APPENDIX A 



MODEL 941 SYLVANIA LASER SPECIFICATIONS 

Output wavelength: 10.6 pm 

Power output : 3 watts 

Beam divergence : 4 mrad 

Beam diameter to 1/e point: 4 mm 

Beam polarization: vertical 

Cavity length: 45 cm 

Bore diameter: " 6 mm 

Cooling requirements: tap water @ 15 to 30 ml/s 

Output mirror: 70% reflecting, dielectric flat 

Piezoelectric mirror: 100% reflecting, 300 cm radius of 

curvature 

Piezoelectric ^ parameter 

As given by Sylvania: -4 x 10 J pm/volt 

Measured: -3.8 x 10 pm/volt 

Laser voltage: 14 kV @ 7.6 ma 

Tube voltage: 9.5 kV @ 7.6 ma 

Series ballast resistance: 600 kfi 

Equivalent tube resistance @ 7.6 ma: 1.27 MQ 

Pressure broadened linewidth: 100 MHz 

Doppler broadened linewidth: 50 MHz 

jj- mode separation: 333 MHz 

Cavity linewidth: 19 MHz 
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APPENDIX B 



STARK CELL SPECIFICATIONS 

Two parallel plates totally encased in plastic with a 
3 mm x 10 mm light beam channel at the axis of the cell 
between the plates. Channel is sealed off with one inch 
NaCl windows at each end and is vacuum tight. Channel is 
filled through one port. 

Plate separation: 7 mm 

Plate width: 3.6 cm 

Plate length: 20.3 cm 
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APPENDIX C 



HP 9810A CALCULATOR PROGRAM 
FOR THE VOIGT PROFILE 

The Voigt profile is the convolution of a Lorentzian 
and a Gaussian profile. This HP 9810A program convolves 
numerically a Gaussian g^ (x) and a Lorentzian g ? (x) which 
have equal heights at the center frequency, i.e., 



MV - MV • (ci) 

Solving (Cl) and setting v Q = 0 and = \ the Gaussian 
and Lorentzian profiles are simplified to 



and 



g 1 (x) = exp [-ttx z ] 



g 2 (x) 



[1 + (TTX) 2 ] 



(C2) 



These were the equations plotted in Fig. 5. Equations (C2) 
are convolved according to 



x 

‘ max 

g( - x ' ) = x / g 1 (y) g 2 U * y) d y 

min 



(C3) 
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The flow diagram and complete program are given on the follow- 
ing pages. The counter containing the value of n is 
register (000). The value of y is stored in(b), the value 
of x in (a). dx and dy are stored in (100) . 1 /tt is stored 

in (099) . (098) is a temporary storage for holding the 
previous value of f for the numerical integration routine. 
Registers (001) through (093) are reserved for storing S n 

where S goes into (n) . x is stored in (097) . 
n max 

_ j na x min i stored in (096) . n + 1 is stored 

max dx 

in (095) and A which is the normalizing factor is stored 
in (094). All registers should be cleared before beginning. 
Begin by prestoring these typical values: 



dx = dy = .1 -»■ (100) 



x = -x . = 4 (097) 

max mm 

n = 80 + (096) 
max ^ 



y = -4 -> (b) 
x = -4 -*■ (a) 



(^o 2 = \ - (099) 

7T 

Set up the plotter for x max = 4, x min = -4 and y max = 1, 

v =0 Load the program and press continue. The computer 
7 min 

will plot out the Voigt profile in approximately 20 minutes. 
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APPENDIX D 



HP 9810A CALCULATOR PROGRAM FOR 
STARK FREQUENCY DISCRIMINANT 

I. FREQUENCY DISCRIMINANT FOR PYROELECTRIC DETECTOR 

The following program plots t> n in Eq. (66) as a function 
of 6, where f(t) is given by (64) with S = 1 and Av = 1. 

The first step in the program, K2, sets the calculator to 
work in radians. 5 is prestored in (000), Av ’ is prestored 
in (006) , and n (which determines the harmonic) is pre- 
stored in (007). The values of the constants A, B, C, D, 
and E in (64) are calculated and stored in registers (001) , 
(002) , (003) , (004) , and (005) , respectively, in the 
program. The variable t is stored in (a) and dt in (b) . 
(008) is a temporary store for the numerator in (64) . 

(012) is a temporary store for the numerical integration 
routine. (010) accumulates the value of the integral in 
( 66 ). 

Begin with all registers cleared and the program loaded. 
Set up plotter for -1 < x < 1 and -2 < y < 2 or some other 
appropriate limits. Prestore -1 in register (000), .5 in 
(006), 1 in (007), .01 or .05 in (b) . Press continue and 
Fig. 20 will be plotted out from left to right. 
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II. FREQUENCY DISCRIMINANT FOR HgCdTe DETECTOR 

The following program plots in Eq. (65) as a function 
of 6, where f(t) is given by (61) with S = 1 and Av = 1 . As 
before (000) , (006) , (007) , and (b) contain the appropriate 
prestored values. In the program (Oil) is the numerical 
integration routine temporary storage register and (009) 
accumulates the value of the integral. 
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APPENDIX E 



LOW PASS L-SECTION FILTER 

Below is a schematic of the simple filter added to the 
high voltage D.C. laser power supply to reduce the 60 Hz 
current ripple. The filter consisted of a .5 microfarad 
16 kilovolt capacitor and a 45 Henry choke which was 
actually the secondary winding of a 13 kV transformer. 

Since the input voltage was 15 kV, the transformer, being 
insulated for 2000 volts less, was isolated from ground on 
plastic supports. The filter successfully reduced the 60 Hz 
ripple in the power supply output by a factor of 50. 
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APPENDIX F 



COPPER CONSTANTAN THERMOCOUPLE TEMPERATURE VS. E.M.F. 



The following table gives values of thermocouple 
temperature and voltage with reference to ice water. 



E.M.F. (mV) 



Temperature (°C) 
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1.9 
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10.28 
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30.11 
32.56 
35.00 
57.39 

39.78 

42.11 
44.50 

46.83 
49.17 
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APPENDIX G 



RELATIONSHIP OF TUBE WALL TEMPERATURE 
TO COOLING WATER OUT TEMPERATURE 

The laser discharge region between the cavity mirrors 
is a pyrex glass tube .7 cm O.D. and 34.3 cm long sur- 
rounded by a glass water jacket 3 cm I.D. It is desired 
to find the average wall temperature of the .7 cm discharge 
tube knowing the bulk water temperature out. Reference 24 
is used for this calculation. The problem is an annular 
tube flow type. It is assumed that the outer water jacket 
tube is at the bulk water temperature and that no heat 
flows through it to the surrounding air. 

The following parameters are given, through experimental 
measurement or a handbook of physical constants. 

flow rate: 14.6 ml/s 

inlet temp . : 22 . 1°C 

outlet temp.: 23.8°C 

voltage across tube and balast resistor: 14 kV 

voltage across tube: 9.5 kV 

tube current: 7.6 mA 

laser power output: 3 watts 

specific heat of water: Cp = 1 cal/g-C° 

density of water: p = 1 g/cm" 5 

viscosity of water at 22°C: u = .009403 g/cm-s 
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thermal conductivity of water at 0°C: k w = .001348 cal/cm-s-°C 

thermal conductivity of water at 40°C: k = .001499 cal/cm-s-°C 

w 

thermal conductivity of pyrex: k g = .0026 cal/cm-s-°C 

hydraulic diameter 2(r Q - r i ): D g = 2 . 3 cm 

tube surface area" A s = tt ( . 7 cm) (34.3 cm) = 75.4 cm 2 

1. Determine if the flow is laminar or turbulent by calcu- 
lating the Reynolds number. 

2 2 2 

cross sectional area of water: A = ir(r - r. ) = 6.68 cm 

o i 

velocity of water: V = 14,6 ml ^ s = 2.19 cm/s 

6.68 cm 

Re - = (1 = 536 

U .009403 g/cm-s 

. . Flow is definitely laminar. 

2. Determine Nusselt number for laminar flow in a circular 
tube annulus, using [24]. 



r . 




Nu i = 8.212 

3. Find the heat added to the water by the laser discharge 
tube only. This is approximately equal to the voltage 
across the tube times the current minus the power output. 

^ =(9.5 x 10 3 ) (7 . 6 x 10"- 3 ) - 3 = 69.2 watts = 16.54 cal/s 
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4 . 



Find the temperature difference AT across the tube and 
the bulk water temperature given the temperature out. 




16.54 cal/ s = (1 |=jU0 AT (14.6 g/s) 



AT = 1.1 C° 



T + T . T + fT 
o l _ o v o 



1 . 1 ) 



T 



1.1 



b 2 



2 



o 



5. The surface coefficient of the inner tube is found 
from the Nusselt number according to 




6. The rate of heat transfer through a . 5 mm glass tube 

wall with an inside surface temperature Tw^ and an out- 
side surface temperature Tw s is given by 




16.54 = (.0026 



cm-s- c 



cal 

rr 



0(75.4 cm 2 )(Tw i - Tw s )/.05 cm 



Tw. - Tw = 4.2 C° 
1 s 
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7. Derive a general formula for the average tube wall 

temperature given that the rate of heat transfer through 
the tube wall equals the rate through the surface 
layer equals the rate dissipated by the water. 
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From the above it is easily shown that 
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Substitution of the equations for T^ and h s in 4 and 5 
into this equation yields 
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Clearly this proves that ATw^ = AT q . 
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APPENDIX H 



PHOTOGRAPHS OF EQUIPMENT 

Figure (H-l) is a photograph of the whole experimental 
setup. The L-section filter is on top of the laser power 
supply in the lower right-hand corner. The cooling 
water temperature control, filter, and regulator are at 
the sink in the upper right-hand corner. 

Figure (H-2) is a photograph of the oscilloscope and 
rack-mounted instruments. From top to bottom the instru- 
ments in the rack are: audio power amplifier, lock-in 

amplifier, 3kV D.C. power supply, 2kV D.C. power supply, 
audio oscillator, and nonovoltmeter . 

Figure (H-3) is a photograph of the laser and optical 
components. The laser is middle right in the photograph. 
Note the thermometer and thermocouple in the coolant tube 
coming out of the rear of the laser. The gas bottles and 
dewar for the thermocouple ice bath are at the bottom of 
the photograph. From the laser counterclockwise are: the 

chopper, beam splitter, spectrum analyzer. Stark cell, and 
pyroelectric detector. 

Figure (H-4) is a photograph of the Stark cell. From 
left to right are: the beam splitter, chopper motor, laser, 

Stark cell, and detector. Note that the chopper motor is 
acoustically isolated from the table with a metal bar sup- 
ported directly to the floor. 
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FIGURE H-l 






FIGURE H-2 
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